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TeV Gamma-Rays from Old Supernova Remnants 
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ABSTRACT 

We study the emission from an old supernova remnant (SNR) with an age of 
around 10^ yrs and that from a giant molecular cloud (GMC) encountered by the 
SNR. When the SNR age is around 10^ yrs, proton acceleration is efficient enough to 
emit TeV 7-rays both at the shock of the SNR and that in the GMC. The maximum 
energy of primarily accelerated electrons is so small that TcV 7-rays and X-rays are 
dominated by hadronic processes, 7r°-decay and synchrotron radiation from secondary 
electrons, respectively. However, if the SNR is older than several 10^ yrs, there are 
few high-energy particles emitting TeV 7-rays because of the energy loss effect and/or 
the wave damping effect occurring at low-velocity isothermal shocks. For old SNRs 
or SNR-GMC interacting systems capable of generating TeV 7-ray emitting particles, 
we calculated the ratio of TcV 7-ray (1-10 TeV) to X-ray (2-10 keV) energy flux and 
found that it can be more than ~ 10^. Such a source showing large flux ratio may be 
a possible origin of recently discovered unidentified TeV sources. 
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1 INTRODUCTION 

The most probable cosmic-ray accelerator in our Galaxy is 
the young supernova remnant (SNR). The detection of syn- 
chrotron X-rays from shells of young SNRs provides us the 



At present, the leptonic process is 
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J.. iKovama et al.lll995l 
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iBamba et alJ 


E003a b. 
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|b|). So far, the evidence for hadron acceler- 



ation, however, has not yet been obtained. High energy 7- 
ray observations may give us important informati on on the 
accel e rated protons l|N aito fc Takaharal [T9 94: Dr urv et alJ 
I1994I: lAharonian et alJ Il994l : lAharonianfc AtovanI 119961) . 
For example, TeV 7-rays are detected from the SNRs 
RX J 1713.7-3946 jEnomoto et alj I200I lAharonian et alJ 
2004a^ and RX J0 852.0-4622 jKatariri et alJ 12005 " 



Aharonian et alJl2005cl) . which can be originated in either 



the decay of neutral pion, arising from the collision of high 
energy protons and interstellar matter, or CMB photons up- 
scattered by accelera,ted electrons (e . g. iPannuti et al]l2003 ; 
iLazendic etall 120041: lEUisonl l200lt iBamba et alJ 120051 : 



lUchivama et al.ll2005^ . 

not yet ruled out. 

Recently, a survey of the inner part of our Gala xy has 
revea l ed sev era l new TeV 7-ray sources fAharonia n et alJ 
I200I l2005allbl I2OO6I) . For some of them, no counterpart 
has been found in any other wave lengths yet. Their spec- 
tra are rather hard, and the flux above 1 TeV is around 
1 X lO"'^^ cm~^s~^. They are extended with the angular 
size of around 0.1°. They should be galactic origin because 
all are located along the galactic plane. In Table we list 
the properties of the Galactic TeV 7-ray sources that are 
also observed in the X-ray bands (together with a typ- 
ical young SNR SN 1006). One of the newly discovered 
sources, HESS J1303— 631 tAharonian e t a l. 2005b') . was ob- 
served by Chandra, and no ob vious counterpart was revealed 
jMukheriee fc Halpernll2005^ . Then, the flux ratio, defined 
by 
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TcV/X 



F^(l - 10 TeV) 
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Fx(2-10keV) ' 

is more than ~ 2. On the other hand, young SNRs, Gas A, 
RX J1713.7-3946 and RX J0852.0-4622, have 7?Tev/x less 
than ~ 2. Although at present the lower limit on -RtcV/x 
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Table 1. Observed properties of Galactic TeV 7-ray sources and a typical young SNR, SN 1006. 



Name 


F^{1-10 TeV)'* 


Fx (2-10 keV)^ 


-RtcV/x'^ 


Classification'^ 


References^ 


Crab 


56 


2.1 X 10" 


2.6 X 10-3 


PWN 


(1), (2) 


ROW 89 


16 


38 


0.41 


PWN 


(3), (4) 


Cas A 


1.9 


1.3 X 10^ 


0.014 


SNR 


(5), (6) 


SN 1006 (NE) 


< 2.1 


19 


< 0.11 


SNR 


(7), (8) 


RX J1713. 7-3946 


35 


5 X 10^ 


0.07 


SNR 


(9), (10) 


RX J0852.0-4622 


69 


> 32 


< 2.1 


SNR 


(11), (12) 


HESS J1813-178 


8.9 


7.0 


1.3 


SNR 


(13), (14) 


HESS J1640-465 


7.1 


0.37 


19 


SNR 


(13), (15) 


TeV J2032+4130 


1.9 


< 2 


> 1 


unID 


(16) 


HESS J1303-631 


10 


< 6 


> 2 


unID 


(17), (18) 



°' Derived for best fitted values (in units of 10 ^^erg s ^cm 
Derived for best fitted values (in units of 10" ^'^erg s" ^cm~ 
(1-10 TeV)/Fx(2-10 keV) . 



Aharonian ct al. 1 2004b): (2) Willi ngale et al.N200j) : f3')!Aliaronia n et al.l i200 5c'): (4) DeLanev et al 



(5') IAliaronian et al.l 120011): ('6'llAllen et al.lJl997D: _^72^3iaronian ct al.'('2005d); (8) Ozaki & Kovam; 



(91 lAharonian et alj ^2004a^: flO l iPannuti et al . (2003); (11) Aharonian ct aj^ 

2001); (13) Aharonian ct al. 12006); (14) Broean ct al. 12005); (15) Sugizaki et alJ 
C17) .Aharonian ct al.. 12005b.) : (18) .Mukheriee fc Halpcrn. I.2005i) . 



(12)|glanc ct al, 
(16) lAharonian et al.. 1.2002.) 



for HESS J1303— 631 is comparable to that of young SNRs, 
it may become much larger with forthcoming deeper X- 
ray observations. Furthermore, HESS J 164 0— 465 has an X- 
ray counterpart that is identified as SNR iAharonian et all 
[2OO6). This object has i?Tov/x ~ 19. Unhke young SNRs, 
these TeV sources with large RtbV/x may show evidence 
for hadron acceleration because inverse-Compton scenario 
requires unusually small magnetic field strength {<C 1 /iG) . 

A large value of RtcV/x is expected for old SNRs. As 
the SNR ages, the shock velocity decreases. In general, pri- 
mary electron acceleration is limited by synchrotron cooling. 
Then, the roll-off energy of electron synchrotron radiation 
is much smaller than that of typical young SNRs, so that 
small synchrotron X-ray flux is expected (e.g.. lSturner et all 
^93). It is also important to consider the association with 
a giant molecular cloud (GMC). Because of large volume, 
old SNRs may encounter the GMC. In this paper, we study 
how large RtcV /x can become for single SNRs and the SNR- 
GMC interacting systems. 



2 EVOLUTION OF SNR 

We consider a simple analytical model of the shock dynamics 
of SNRs expanding into the uniform ambient medium with 
the density no. The SNRs evolve through three phases: the 
free expansion phase, the Sedov- Taylor phase, and the ra- 
diative phase. We assume the shock velocity Va as a function 
of the SNR age tagc as 

{Vi (0 < tagc < tl) 

Vi{Uge/tl)-^^'^ (tl < tagc < t2) , (2) 

«i(t2Al)"^/^(tagc/t2)"^/^ {t2 < tagc) 

where ti = (3£/27rmHnowf )^/^ = 2.1 x 
lO^iEii/noy/y-^^^ yrs and ta = 4 x lO^'Sg^^n"'"'/" yrs 
jBlondin et al.l ^9^, and Vi = Vi^r,lO^ cm s~^ and 



E = £51 10 ergs are the initial velocity, and the initial 
energy of the ejecta, respectively. Here we adopt the 
expansion law of the radiative phase, Vs oc fZ^J^ , instead of 

— 5/7 

Va OC tago , bccause the former gives a b etter approximation 
in the epoch around tap e ~ 10^ yrs (iBlondin et al.l [19981 : 
iBandiera fc Petrukl20o3) . The evolution of the shock radius 
Ra — J Vadt is also calculated. In this paper, we adopt 
E51 = «i,9 = no = 1 as a fiducial case. Then, we find 
Ra = 91 pc and Va = 1.1 x 10^ cm s~^ at tage = 3 x 10^ yrs. 
Note that the effect of the magnetic field on SNR dy- 
namics can be neglected unt il fagc of several 10^ yrs (e.g., 
iHanavama fc Tomisakall2006^ . 



3 EMISSION FROM AN SNR 

At first, we consider the spectrum of emitting high-energy 
particl es that get their energy via diffusive sh ock acceler- 
ation JPrurvl Il983l: iBlandford fc Eichleil Il987*) . Generally 
speaking, in order to obtain the energy spectrum of accel- 
erated particles, time-dependent kinetic equ ation should be 
solve d as investigated by many authors (e.g.. lBerezhko et alJ 
Il996l) . eve n including the nonline ar effects via accelerated 
particles llMalkov fc Drurvl 1200]]) . Such calculations well 
match the observational facts when free parameters such as 
the injection rate and the magnetic field are chosen appro- 
priately. Instead, we assume, in this paper, that at an arbi- 
trary epoch, the spectral form is a power-law with the index 
p and the exponential cut-off at Emax, i.e., oc E~^e~^^^'^'^^ , 
where E'max evolving with time. This approximation may 
be valid and be useful to extract the basic properties of 
high-energy emission from a SNR because of the following 
reasons. High-energy particles are produced by the diffu- 
sive shock acceleration and suffer adiabatic expansion af- 
ter they are transported downstream of the shock and lose 
their energy. Hence, at the given epoch, the spectrum for 
high-energy particles is dominated by those which are being 
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accelerated at that time, in other words, the energy spec- 
trum of particles does not so much depend on the past ac- 
celeration history. Especially, we are now interested in the 
energy region near the upper end of the spectrum because 
as seen in the following, both X-rays and TeV 7-rays are 
produced by particles with energy near i?max- In this energy 
regime, our assumed form of the spectrum may be a good 
approximation . 

Once the SNR dynamics is given, the maximum energy 
of accelerated particles is calculated. The maximum energy 
of accelerated protons, -Emax.p, is determined by 



The wide-band radiation spectrum 



min{tagc,tpp} 



while that of accelerated electrons, 

^acc — min{tagc , ^synch } , 



(3) 

is determined by 
(4) 



where tage, tacc, tsynch, and tpp are the age of the SNR, the 
acceleration time scale, the synchrotron loss time scale, and 
the pion-production loss time scale, respectively. Since we 
assume that the energy spectrum does not depend on the 
past history, it is also a good approximation to determine 
i?max using Eqs. or @. Indeed, at least for young SNRs, 
one can confir m our estimation of -E max coincides with simu- 
lation results teerezhko et alj |2002ll . For the diffusive shock 
acceleration, the acceleration time is given as 



_ 20fec£n,ax 
t-acc 7-1 ■> ' 

eBdvi 

where 

_ 0.05r(/ + rg) 



(5) 



(6) 



and r is the compression ratio, and / and g are functions of 
the shock angle 9 and gyro-factors jyu and r/d that are given 



f{r]u,0)=r]u{cos + r sin 



Nl/2 



g{i]d,d) = r]d{cos^ 6 + sin' 6) ^ 



, (7) 



(8) 



respectively jYamazaki et al .l2004l : IJokipill987l) . The down- 
stream magnetic field is given by = rSisM, where 
BiSM ~ 10-BisM,-5 fJ,G is the ISM magnetic field and 
we adopt Bism.-s ~ 1 as a typical value. For high- 
energy protons, the energy loss time scale through pion- 
production is given by tpp = 5.3 x 10^n~^ yrs, where n 
is the number density of the acceleration site. For elec- 
trons, the synchrotron loss time scale is given by fsynch = 
1.25 X 10''(i;max,e/10TeV)-^(Sd/10/iG)-2 yrs. As long as 
tage < tpp, the proton acceleration is age-limited, and we 
obtain 



Smax.D = 1.6 X 10^ TeV h-\', ' 



fa 



(9) 



10 /iG y V lO^yrs 

where Vs,8 = Us/lO^cm s^^. The electron acceleration is loss- 
limited when fagc 10'^ yrs, and then we derive 

Bmax.e = 14 TeV h-^^'vs.H i I 



' yiOpG J 



(10) 



The minimum energy, i5min,c and i5min,p are simply taken as 
the rest-mass energy of electrons and protons, respectively. 



given En 



and En 



at a certain fa 



is calculated for 
. We consider ra- 



diation processes from primary electrons: the synchrotron, 
inverse- Compton (IC), and bremsstrahlung emissions, and 
from primary protons: tv^ decay 7-rays, and synchrotron and 
bremsstrahlung emissions from secondary electrons arising 
from the decay of charged pions. We assume the electron- 
proton ratio, that is the ratio of the electron distribution 
function to that of protons for a fixed energy in relativistic 
regimes, of Kcp = 1 x 10~^. For average cosmic-rays in our 
Galaxy, Kcp is about 10~^, which is an order of magnitude 
larger than our adopted value. However, in general the aver- 
age value may be different from that at the acceleration site 
because of the propagation effect. Indeed, as seen in the fol- 
lowing, if we adopt K^p = 1 x 10~^, theoretically calculated 
flux ratio -RtcV /x for a young SNR is similar to the observed 
one. For old SNRs, which we are most interested in, the re- 
sult for -RxeV/x is unchanged even if K^p varies more than 
one order of magnitude because accelerated electrons do not 
contribute to TeV and X-ray emission. 

Primarily accelerated electrons produce synchrotron 
emission with the roll-off frequency given by v^qu ~ 1.6 x 



Bd/10MG)(Sn,ax,c/10TeV)^ Hz jRevnolds fc Keohan3 
Using Eq. IIOII . i^roii can be rewritten as 

iv.oii ~ 3.2 X lO^'^h'^vls Hz , (11) 

in the loss- limited case (see also lAharonian fc Atovarll999t) . 

Note that UroU does not depend on Bism- 



3.1 Emission from a young SNR in the Sedov 
phase 

As can be seen from Eqs. @ and llOII . once the SNR dy- 
namics and i3isM(= B^/r) are fixed, only necessary are the 
values of r and h in order to calculate -Emax.p and i?max,o for 
a given time. In the Sedov phase (fi = 2 x 10^ yrs < fagc < 
f2 = 4 X 10* yrs for £51 = no = Vi^Q — 1), the shock of an 
SNR is strong and adiabatic, so that the density compression 
ratio, r, is near 4. Then in the Bohm limit case, f^u ~ f?d ~ 1, 
we find /i ~ 1 for arbitrary 9. When the nonlinear ofi^ccts are 
considered, r becomes as large as ~ 7 (|Bcrczhko ct al. 200^, 
however, in this case, one can find that h weakly depends on 
r. Hence in this paper, we neglect the nonlinear effects for 
simplicity. When h is around unity, we can reproduce the 
observed value of Vy-ow for young SNRs (e.g., iBamba et alJ 



2003b, 2005a b). Hence, for fago < 10" yrs, we adopt r = 4 
and /i = 1 as a fiducial value. The power-law index of ac- 
celerated particles is fixed as p = 2.2, which is typical for 
young SNRs. 

Here we assume that an SNR stores energy, 10^° ergs, of 
high-energy protons. As long as accelerated particles main- 
tains upstream turbulence via streaming instability, the self- 
confinement of accelerated particles is efficient, so that the 
escape of them upstream may be neglected. Furthermore, ra- 
diative efficiency is also small. As discussed previously, the 
energy of accelerated particles is reconverted to the expan- 
sion energy through the adiabatic expansion, which implies 
the expansion energy remains almost unchanged. Therefore, 



^ Seell http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec| 

for the erratum of coefBcients. 
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Figure 1. vFi, spectrum of a single SNR with an age 1 X 10^ yrs 
that stores energy, 10^" ergs, of high-energy protons. The thick- 
solid line shows the total non-thermal flux. Hadronic emis- 
sions are Tr'^-decay 7-rays (thin-solid), synchrotron (dot-dashed) 
and bremsstrahlung emission (short-and-long dashed) from sec- 
ondary electrons produced by charged pion. Leptonic emissions 
are synchrotron (long-dashed), inverse-Compton (dotted) and 
bremsstrahlung (short-dashed) emission by primary electrons. 
The source is located at 1 kpc. We adopt £51 = tii^g = no = 



h — -BlSM,-5 
for details). 



1, p = 2.2, r = 4 and 




1 X 10"'^ (see text 



10 15 



10 -5 5 10 15 

LogjE/eM) 



Figure 2. The same as Fig.^but for an SNR age of 1 X Iff* yrs. 



Meanings of each lines are the same as those in Fig. The ab- 
solute value of the observed flux is uncertain, hence the flux is 
arbitrary scaled. We adopt £51 = i)i g = riQ = h = Bism —5 = li 



1.5, and Kc 



1 X 10 (see text for details). 



if the injection rate is time-independent, the energy of high- 
energy protons is also time-independent. 

When tage = 1 X 10^ yrs (then, ti < tagc < ^2), we find 
-Bniax.p = 96 TeV, E^^^,, = 27 TeV and f.oii = 4.8 x lO'^ Hz 
for the fiducial parameters, and find that nonthermal X-rays 
are dominated by the synchrotron emission from primary 
electrons (see Fig.0. TeV 7-rays are dominated by Tr'^-decay 
process. The flux ratio is RtcV/x = 7.6 x 10~^, which is 
consistent with observations for young SNRs 

When tage 
Sedov phase (ti < 



1 x 10* yrs (Fig. 13, an SNR is still in the 
age < t2). For the fiducial parameters. 



En 



61 TeV and E, 



max,c 



6.9 TeV are derived. Also in 



this case, the X-ray band is dominated by the primary syn- 
chrotron radiation. As the SNR ages, froii becomes small, so 
that synchrotron radiation fiux in the X-ray band becomes 
small. Hence the flux ratio becomes large, -RtoV/x = 1-6. 



3.2 Emission from an old SNR in the radiative 
phase 

When an SNR enters into the radiative phase (tagc > ^2), the 
cooling effect becomes important, and the accumulated gas 
forms a dense, cool shell with the number density rish = rno- 
Rewriting Eq. Q as 



Vs = 2.3 X WE. 



ta 



-2/3 



(12) 



lO^yrs 

we find Vg ~ lO'^ cm s-i when tago ~ 10^ yrs. Note that Vg 
is independent of Vi. 

In the radiative phase, the shock is isothermal. Hence 
there is a concern that since the downstream temperature is 
low, the neutral component appears, and Alfven wave turbu- 
lence that scatters nonthermal particles is significantly sup- 
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pressed via the ne utral-ion friction, makinK the accelera tion 
process inefficient dOrurv et al .llQQdiBvkov et all200d) . On 

the other hand, when the gas remains fully ionized around 
the shock front, diffusive shock accelera tion works well in the 
radiative phase. IstiuU fc McKed l|l979^ found that as long as 
Vs,7 > 1.1, where Vs,7 = Vs/lO'cva s~^, UV flux radiated in 
the downstream region is sufficient enough to keep upstream 
state fully ionized. For no = iSsi = 1, we find «s,7 = 1.1 at 
^age — 3 X 10 yrs. Therefore, when fago < 3 x 10 yrs, diffu- 
sive shock acceleration works well, and we can use Eqs. ® 
and lO in order to estimate the maximum energy of accel- 
erated particles. 

The estimation of downstream magnetic field, Bd, 
which appears in Eqs. @ and llOII . is different from that 
for the young SNRs. For no = -Bism,-5 = 1, the magnetic 
pressure is larger than the gas pressure in the upstream re- 
gion. Then, solving the shock jump conditions for the quasi- 
perpendicular [6 ~ 90°) isothermal shock, the compression 
ratio is given by 

r ~ ^2-^ = 6.5 BisM,-5^^nJ^Vs,7 , (13) 

VAl 

where v ai = -BisM(4 7rmj:f no)~^^^ is the upstream Alfven 
velocity l)Spitzeill97Sh . Therefore, the downstream magnetic 
field Bd = rBisM = {SnrnHnoy^^Vs is calculated as 

Bd = 65 /iG nl^'^vs,? ■ (14) 

Recently, using two-dimensional MIfD simulation, 
[fjanavama & Tomisaka (2006) have shown that in the 
case of no — 0.2 (or 1) and £51 = Bism = 0.5, the compres- 
sion ratio at the quasi-perpendicular shock region is about 
7 at tagc of several 10''' yrs, which is roughly consistent 
with our estimation. Using Eqs. (|5J, llOH . 11211 and I14II . we 
obtain 

Smax.p = 34 TeV h''ElY'\-''^"v%^ , (15) 
S,„ax,c = 0.55 TeV ft-^'^n-'/X^^' . (16) 

Note that Eq. 1141 1 is derived for the limiting case of a small 
upstream gas pressure compared to the magnetic pressure 
there. If the explosion occurs where the gas pressure is com- 
parable to or larger than the magnetic pressure, the effect 
of the gas pressure should be taken into account. Then, 
the compression ratio becomes large iSpitzeJ Il978f) . Fur- 
thermore, when the accelerated protons are stored around 
the shock front, a part of their energy goes into the mag- 
netic field energy, c ausing the magnetic field amplification 
jLucek fc Bellll200(i) . Therefore, Eq. lO gives, in fact, the 
lower bound of B^, so that, Eqs. 1151 and 11611 give lower and 
upper bounds of Smax,p(oc Bd) and Bmax,e(oc Bd~^^^), re- 
spectively. In the following, however, neglecting such effects 
for simplicity, we use Eqs. 1141 . 1151 and 1161 to determine 
Bd, Bmax.p and -Emax.c. Oue should also remark that Eq. 1141 
is derived for the quasi-perpendicular shock {6 ~ 90° ) . Then, 
one can consider the accelera t ion at the quasi-perpendicular 
shock as discussed in lJokipiil I^S^)- However, even if ^ 
upstream, then the accelerated particles penetrate into up- 
stream region and generate waves via streaming instability, 
resulting in the generation of oblique or quasi-perpendicular 
magnetic field upstream, because the magnetic field compo- 
nent of generated waves is generally perpendicular to the 
shock normal and can be as large as the original upstream 
magnetic field. 



When tagc = 3 X 10^ yrs, we find r = Bd/BisM ~ 7. 
Hence, in the case of Bohm limit ?7u ~ 77d ~ 1, we find h ^ 1 
as well as for the young SNRs. Until this epoch, the particle 
acceleration proceeds, so that the upstream turbulence is 
maintained and the Bohm limit diffusion can be expected. 
Therefore, we again take /i = 1 as a fiducial value. In the test 
particle approximation, the power-law index of accelerated 
particles p and the compression ratio r are related as p = 
{r + 2)/{r-l) jBlandford fc Eichler.,.,198l . If we take r ~ 7 
as a typical value, p becomes about 1.5. Hence, we adopt 
p — 1.5 as a fiducial value. However, as can be seen in the 
following, our conclusion does not depend on p. 

The total energy of accelerated protons stored in the 
SNR is somewhat uncertain. If the confinement of acceler- 
ated particles works well until t^gc ~ 10^ yrs, the energy 
of accelerated protons may not be so small compared with 
the young SNR case. Otherwise, it is much smaller than 
10^" ergs. In this paper, since we are mainly interested in 
the value of RtbY/x, the total energy of accelerated parti- 
cles, which only determines the normalization of the radiated 
spectra, is not essential for our arguments. Further discus- 
sion on determining the amount of accelerated particles will 
be seen in § |S] 

Figure|3shows the spectrum for tage = 3 x 10^ yrs with 
fiducial parameters (i.e., Bd = 72 /xG, Bmax.p = 42 TeV and 
£max,o = 0.58 TeV). Then the fiux ratio is RtcV/x = 82. 
Since -Emax.o is small, TeV 7-rays via IC and bremsstrahlung 
are suppressed, and therefore, the TeV 7-rays come from 
the 7r°-decay. The roll-off frequency of synchrotron radiation 
from primary electrons is so small, VroU ~ 4.0x10^"' Hz, that 
the secondary synchrotron radiation dominates the X-ray 
band. Also in this subsection, we have assumed the electron- 
to-proton ratio as K^p — Ix 10^^. However, even if we adopt 
much larger K^p, RtcV/x remains unchanged because as can 
be seen in Fig. |3 both X-ray and TeV 7-ray bands lie above 
the cutoff of emission components originating in primary 
electrons. 

We discuss how RtcV/x varies for different values of 
parameters. At first, changes of Bism and Vi do not af- 
fect the conclusion because Eqs. 1141 . 1151 and 1161 tells 
us that Bd, -Emax,p and i5max,c, which determine the overall 
shape of the emission spectrum, do not depend on them. As 
stated above, uncertainty of -K'ep does not affect the value 
of i?TcV/x- We also find that -RtoV/x only weakly depends 
on h and E51 for reasonable parameter ranges. Hence we 
discuss the dependence of no and p in the following. Let 
us first discuss the case of no = 0.5 with other parame- 
ters unchanged. We consider the epoch in which Vs,7 ~ 1.1 
(tage ~ 3.8 X 10^ yrs for our parameters) when -RxeV/x be- 
comes maximum since it increases with time. Then, we find 
Bd = 51 ^G, £„iax,p = 39 TeV and _En,ax,c = 0.70 TeV, re- 
spectively. Again both the X-rays and TeV 7-rays are domi- 
nated by the hadronic emission. Compared with the fiducial 
case, the secondary synchrotron X-ray emission is dim be- 
cause Bmax,p and Bd are small (see the dashed line in Fig.|lJ, 
so that the flux ratio becomes large -RxeV/x = 1-3 x 10^. 
Next, we consider the case p = 2 with other parameters be- 
ing fiducial (see the dotted line in Fig. |1J . TeV 7-rays are 
produced by protons with energies of ~ 10 TeV, while X-rays 
come from secondary electrons with energies of ~ 30 TeV 
originated in ~ 100 TeV protons. Then, for a large p, 
there are less protons with ~ 100 TeV and the dimmer 
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Figure 4. uF,^ spectrum of a single SNR with an age 3 X 10^ yrs 
for various parameter sets. Solid line is the total flux for the 
fiducial set of parameters (ng = 1, p = 1.5, = 72 fiG, 
Emax.p = 42 TeV and i?max,c = 0.58 TeV), and is the same as the 
thick solid line in Fig.|21 Dashed and dotted lines are for no = 0.5 
(Ba = 51 fj.G, Smax.p = 39 TeV and Emax.e = 0.70 TeV) and 
p = 2 with the other parameters fiducial, respectively. 



Figure 5. uFi, spectrum for emission from a shocked GMC that 
is collided by a SNR with an age of 4.6 X 10** yrs. The GMC has 
mass of 10^ Mq, and the meanings of each line are the same as 
those in Fig. ^ The flux is arbitrarily scaled. We adopt E^i = 
Vi,9 = no = h = iJisM,-5 = M5 = 1, Kcp = 1 X 10"^ and p = 1 
(see text for details). 



X-rays. Hence we find the flux ratio becomes as large as 
RtbV /X = 1-6 X 10^. So one can see that when parameters 
are changed within a reasonable range, the value of i?TcV /x 
remains unchanged within a factor of two or three. 

If tago ^ 3 X 10^ yrs, i.e., Vs,7 < 1.1, t hen, the ioniza- 
tion around the shock front is incomplete jShuU fc McKed 
[^T^) and upstream ion-neutral Alfven wave da mping places 
significant restriction o n shock acceleration jPrurv et alJ 
Il996l: iBvkov et alJl200ol) . Once the shock acceleration be- 
comes inefficient, there are few high-energy protons emit- 
ting TeV 7-rays around the SNR shell, because they es- 
cape the SNR shell due to the diffusion; assuming the 
Bohm difi^usion, the escape time for a particle with an en- 
ergy Ed = 10i?cr,ioTcV TeV is estimated as teac ~ 1 x 
10^77"^ (Bd/70/iG)£;^";\oTov^3pc^ yrs, where 77 and A = 



^cr.lOToV^Spc 

SAspc pc are the gyro factor and the thickness of the shell, 



respectively. Therefore, when fa 
are significantly suppressed. 



> 3 X lO^yrs, TeV 7-rays 



o 



10 15 




10 -5 5 10 15 

LogjE/eW) 



4 EMISSION FROM AN OLD SNR 
INTERACTING WITH A GMC 

One may expect that if an old SNR interacts with a GMC, 
TeV 7-ray flux becomes large because of the large density 
in the GMC. Here we consider two cases for the emission 
from old SNR-GMC interacting systems; one comes from a 
shock running into the GMC, and the other from the GMC 
illuminated by particles accelerated at the SNR shock. We 
consider the GMC with the mass Mc = 10^ Ms Mq, and 
the radius Rc = 18 pc Mg''^ iBlitz fc Rosolowskvl liooil . 
Assuming the spherical symmetry, the mean number density 



Figure 6. uF,, emission spectrum from a GMC illuminated by 
high-energy particles arising from the SNR with an age of 3 X 
10^ yrs. The GMC has mass of IO^Mq, and the meanings of 
each line are the same as those in Fig. Q The flux is arbitrarily 
scaled. We adopt E^i = Vi^g = no = h = -Bism,— 5 = A/5 = 1, 
Kep = 1 X 10"*^ and p = 1.5 (see text for details). 



Te V gamma-rays from old SNRs 7 



of the GMC is 



1.7 X lO^Mj ^''^ cm"^. We also assume 

1/2 



that the magne tic field of the GMC is Be = 1 mG n, 
JCrutcheilll99lh . Then we find 

Be = 13 AtG M;^'^ , (17) 

and the Alfven velocity in the GMC is va = 2.2 x 10^ cm s"^. 
In this section, numerical values are for the fiducial param- 
eter set (£51 = «i,9 — no — h — Bism.-s ~ M5 = 1) unless 
otherwise stated. 



4.1 Emission from a shock running into a GMC 

When an SNR shell collides with a GMC with the number 
density ric, the shock front is formed in the GMC. In gen- 
eral, since the geometry of SNR shell and the GMC at the 
collision is unknown, it is hard to precisely determine the 
dynamics of the shell. Here we perform a very rough cal- 
culation. Assuming the momentum conservation, the shell 
velocity, Vs,c, in the GMC can be related with th e velocity, 
Vs, o f the SNR shell just before the collision as JChevalieil 
Il999^ 

(18) 



1 + (ne/nsh)i/2 



The radiative cooling is efficient and the shock is isothermal. 
We find that the velocity of the shock in the GMC is Vs,c ~ 
1.1 X 10^ cm s^^ when the SNR shell with an age of 4.6 x 
10* yrs {vs — 3.9 x 10^ cm s~^ and rish ~ 25 cm~^) collides 
with the GMC with the mass of 10^ Mq. Then, the Alfven 
Mach number of the shell running in the GMC is ~ 50, and 
the compression ratio is estimated as r ^ \/2vs.c/va ~ 70, s o 
that the shock acceleration works well (iBvkov et al.ll200Cl) . 
The downstream magnetic field is 



Bd = (87rmHne)^''^«s,c = 9.1 x 10^ /iC 



(19) 



and we find E^^^^^ = 79 TeV and £max,o = 0.16 TeV. Here 
we adopt p = 1 as a limiting case, because the compression 
ratio is much larger than unity. The simulated spectrum 
is shown in Fig. |S] We assume the electron-to-proton ratio 
Kep = 1 X 10~^ as in § however, our conclusion does not 
change if the ratio becomes ten times larger. Since Bmax.c is 
small, both X-ray and TeV 7-rays are hadronic origin. Be- 
cause of the large magnetic field, the secondary synchrotron 
radiation is strong in the X-ray band, and the flux ratio is 



not so large, R 



ToV/X 



1.9. 



We calculate -RtcV/x for several cases in which the 
shocked GMC has Vs,c = 1.1 x 10^ cm s~^, and find that 
the large magnetic field at the emitting region causes bright 
secondary synchrotron radiation, and that the fiux ratio for 
the emission from shocked GMC does not exceed ~ 20. We 
also note that if we adopt p = 2.0 instead of p = 1.0, with 
other parameters fiducial, we obtain -RxeV/x = 15, which 
implies the value of p does not affect our conclusion. 

If the SNR shell with tagc > 5 x 10** yrs coUides with the 
GMC with Afs = 1, «s,c is smaller than 1.1 x lO'' cm s"\ 
so that the TeV emission from the shock in the GMC is 
significantly weak because of the wave damping effect and 
the energy loss effect as discussed in § 13.21 At this time, the 
density is so high that high-energy protons lose their energy 
via pion production. Hence when SNR with tage ^ 5x 10"* yrs 
collides with GMC, the arising shock in the GMC does not 
emit TeV 7-rays. 



4.2 Emission from a GMC Illuminated by protons 
accelerated at an Old SNR shock 

If an SNR with tagc > 5 x lO"* yrs interacts with a GMC 
with A/5 = 1, particle acceleration at the GMC shock is 
inefficient to produce TeV 7-ray emitting particles because 
Ws,c is smaller than 1.1 x 10^ cm s~^. Nevertheless, the high 
density of GMC works as a target of the high-energy protons 
to produce pions. They are accelerated at the shock front of 
the SNR and penetrate into the GMC, radiating photons. 
We have seen in § EH that the SNR shock itseff can accel- 
erate TeV-7-ray-emitting particles until tage ~ 3 x 10^ yrs. 
Here we consider the emission from the GMC encountered 
by the SNR with tago = 3 x 10"" yrs. The simulated spec- 



= 42 TeV and 



trum is shown in Fig. |S| We use Bir 
Bmax.c ~ 0.58 TeV as seen for the fiducial case in § 13.21 
while we adopt the magnetic field of the emission region as 
Be = 13 /iG. Both the X-rays and TeV 7-rays are again 
hadronic origin. Compared with Fig. |51 the secondary syn- 
chrotron radiation is weaker in the X-ray band because of 
smaller magnetic field. We find BtoV/x ~ 5.9 x 10^. If the 
GMC is encountered by the SNR with tagc > 3 x 10^ yrs, 
particles emitting TeV 7-rays do not exist around the shock 
front as seen in S 13.21 

We calculate Bxev /x for several cases in which parame- 
ters are changed within a reasonable range, and find that the 
value of i?TcV /x remains unchanged within a factor of three. 
We also note that if we adopt p = 2.0 instead of p = 1.5, with 
other parameters fiducial, we obtain RtcV/x = 1.4 x 10'^, 
which implies the value of p does not affect our conclusion. 

For Be = 13 /iG, both TeV 7-rays and X-rays are origi- 
nated in the accelerated protons with energies of more than 
10 TeV. Let S = Spc pc be the separation between the GMC 
and the SNR shell producing high-energy particles. The time 
for a particle with an energy Ecr = lOBcr.ioTcV TeV to 
reach the GMC is estimated by the diffusion time, tdif ~ 
2 X 10^r?-iBisM,-5B- 

r^oTeV^pc^ y^s, whcrc Tj IS the gyro 
factor. If 5pc ~ 1, we expect tdn is much smaller than the 
age of the SNR. Hence all particles in the energy regime in 
which we are interested reach the target GMC almost si- 
multaneously, which implies that the particle spectra at the 
GMC does not so much depend on the character of propa- 
gations from the SNR sheU to the GMC. If Spc > 1, effects 
of energy-dependent diffusion should be considered. Such a 
detailed calculation is not considered here. 



5 DISCUSSIONS 

For an old SNR in the radiative phase, the maximum energy 
of primary electrons is so small that emissions via leptonic 
processes are vanishingly small both in the X-ray and TeV 
7-ray bands. On the other hand, there still exists acceler- 
ated protons with energies of more than ~ 10 TeV. These 
particles emit the TeV 7-rays via 7r°-decay and synchrotron 
X-rays from secondary electrons generated by charged pions. 
We find that the ratio 



7?' 



TcV/X 



10 TeV) 



Fx(2- 10 keV) 



could be more than ~ 10 . Such sources may be an origin of 
recently discovered unidentified TeV sources, and give us the 
evidence for hadron acceleration. We might have to consider 
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Table 2. Calculated flux ratios for various bands. 







-RxeV/x" 


-RtgV/Gov'' 


^TeV/radio'^ 


Fig. 


1 


7.6 X 10-2 


0.36 


9.8 


Fig. 


2 


1.6 


0.25 


6.6 


Fig. 


3 


82 


6.6 


47 


Fig. 


5 


6.9 


74 


9.7 X 10^ 


Fig. 


6 


5.9 X 10^ 


5.8 


3.5 X 102 



1 X 10 "■'nsh,o.86-Ep,5orfkp^c cm ^s" 



lated as F{> TeV) 

where rish ~ lO^'^^nsh.o.se cm-^ and d — dkpc kpc are the 
density of the SNR shell and the distance to the source, re- 
spectively. The maximum distance to the source is dmax ~ 



1/2 



1/2 

h 86^p 50 '^P^- Then the volume fraction of the survey 



region is 

(7r/6)d„ 



7r(10 kpc) 



0.30 nsh,o.86-E'p,5o 



RtoV/x = ^7(1-10 TeV)/Fx{2-10 keV) . 
flTeV/GcV = F^(l-10 TeV)/FGev(l-10 GeV) . 
Rxev/radio = F-i{l-10 TeV)/F,,dio (lO^-lO" Hz) 



the interaction between the GMC and the SNR with an age 
of ~ 10^ yrs. For tagc ^ 10^ yrs, SNR radius is so small that 
there is only a few SNRs interacting with a GMC. On the 
other hand, if the SNR with tagc ^ 10^ yrs encounters the 
GMC, there are few high-energy particles emitting TeV 7- 
rays around shocks of the SNR and the GMC because of the 
energy loss effect and/or the wave damping effect occurring 
at low-velocity isothermal shocks. 

Actually detected TeV sources have the energy 



flux F^{\ - 10 TeV ) 



10"' -10"''erg cm 
Hence if -Rtcv/x 10 



s 

,2-3 



jAharonian et alj|2005al . 12006 

then Fx{2 - 10 keV) > 10"^* erg cm-^s-\ Such diffuse, 
extended source can be detected with current X-ray tele- 
scopes {Suzaku, Chandra, and XMM-Newton). Especially, X- 
ray Imaging Spectrometer (XIS) onboard Suzaku is capable 
of observing dim, diffuse X-ray sources with low background 
in the hard X-ray band. Less energetic protons emit GeV 7- 
rays that might have been detected by EGRET or may be 
detected by GLAST in the future. However, for isothermal 
shocks, the spectrum of accelerated particles may be hard 
(p < 2). Indeed, som e of the newly discovered TeV sources 
show hard spectrum jAharonian et al]l2005al . l2006l) . which 
may imply p < 2. In such a case, GeV emission becomes 
dim as shown in Table |5| where we calculate i?TcV/GeV for 
the case of Figures 1-3, 5 and 6. Radio emission is also ex- 
pected. However, as can be seen in Table|21 the sources with 
large iZxeV/x have large iZxeV/radio, so that the radio emis- 
sion may be dim. This fact might account for a large scatter 
of iixev/radio fouud bv iHelfaud et al? (2005). Furthermore, 
the source is confined in a galactic plane region and radio 
emission may be obscured by diffuse components. 

In this paper, we consider TeV 7-ray emission from (1) 
a single old SNR with iage ^3x10^^ yrs (§Ell, (2) a shocked 
GMC coUided by the old SNR with tago ~ 5 x 10* yrs 
(§ 14. in . and (3) an unshocked GMC illuminated by high- 
energy particles arising at the shock of the old SNR with 
tage ~ 3 X 10^ yrs (i; l4.2L respectively. As in the foUowings, 
we can roughly estimate, for each case, the expected num- 
ber of observed TeV sources, NtbV, which have observed 
flux larger than 3% of the Crab flux (6 x 10"" cm'^s-^ 
above 1 TeV) and lie in the inner part of the galactic plane 
with -30° < Z < 30° fcf. lAharonian et al.l2005al.l2006r) . The 
number of detected TeV sources gives us important informa- 
tion on the energetics of old SNR emission. 

Let us first consider the case (1). Let Ep = 
Ep^soW^^ ergs be the energy of accelerated protons stored 
in the SNR. Then, the observed flux of TeV 7-rays is calcu- 



Hence, for assumed SN explosion rate of 1 x 10 ~r_2 yr 
we obtain 

iVxcV ~ 3 x 10^ r_2^i 

~ 9 X 10^ r_2nsh,o.86-Ep,50 ■ 

Next, we calculate A'xcV for case (3). In order to esti- 
mate the number of GMCs in our Galaxy, we adopt the mass 
function of GMCs, f{Mc) oc M^^'^ with the maximum and 
the minimum mass, Mmax = 2x lO^M^ and Mmin = 30Mq , 
respectively jBlitz fc Rosolowskvl |20(M1 . If the total mass 
of GMCs in our Galaxy is 2 x 10^ Mq, then the number 
of GMCs with mass of around Mc = W'^MsMq is about 
7 X 10^ M~'^^'^. The SNR with an age of 3 x 10^ yrs has a 
radius Rs ~ 90 pc, so that the geometrical factor, represent- 
ing the fraction of accelerated protons colliding the GMC 
and emitting 7-rays, is estimated as a = (Rc/Rs)^ /i ~ 
1 X 10"^. The observed flux of TeV 7-rays is calculated as 
F{> TeV) ~ 3 X 10""a- 



2nc,2.2-Bp,5orfkpc cm s , where 
a — a_2l0"'^, and ric = 10^'^nc,2.2 cm"^ is the den- 
sity of the GMC. The maximum distance to the source is 

1/2 1/2 1/2 

dmax ~ 6a_2 n/y o.Erlm kpc. Then the volume fraction of 



-"-2 "'c,2.2-^p,50 

the survey region is 



/3i 



(7r/6)d„ 



7r(10 kpc) 



7 X 10" 



a-2nc,2.2-C/p,50 



Assuming the volume of the galactic disk of 4 x 10^" pc'' 
the mean separation of GMCs is £ ~ 4 x 10^ pc. Then 
the probability that a SNR collides with a GMC is 



/32 



47r /^\ 
3 W J 



5 X 10"^£400"^ 



Rs 



where 

A^'xov 



90 pc 

400 ^400 pc. Therefore, we derive 



3 X 10'^ r_2/3i/32 

10 r_2Q-2nc,2.2£'p,50^400~ 



Rs 



90 pc 



Hence £"^,50 > 0.1 is required so that A'^xcV ^ 1- 

Finally, we consider the case (2). We assume that the 
energy of accelerated protons is 1 x 10'**i5p 4g ergs, which 
is much smaller than the thermal energy stored in the 
shocked GMC. Then, the observed TeV 7-ray flux is cal- 
culated as F{> TeV) ~ 3 x 10"^nc,d,4.i-Bp_48d^p2^ cm"2s"\ 

where nc,d, = 10'''^nc,d,4.i cm"^ is the density of shocked 
GMC. The maximum distance to the source is calculated as 



66n 



1/2 



1/2 



kpc, which is larger than the size 



of the Galactic disk. Therefore, the volume fraction is about 
/3i ~ 0.3. Since the SNR with an age of 5 x 10* yrs has a 
radius Rs ~ 40 pc, the collision probability is 



/32 



47r 

T 



4 X 10"^^400"^ 



Rs 



40 pc 
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Hence we obtain 



iVTcv ~ 5 X 10^ r_2/3i/32 

~ 0.6 r_2(/3i/0.3)^40o" 



Rs 



40 pc 



Taking into account for the uncertainties for Rs and/or ^400, 
we can expect A'tcV > 1 for reasonable values of parameters. 

At present, the number of unidentified TeV sources 
showing large -RtcV/x is only a few. If the number of ac- 
tually detected TeV sources does not increase (i.e., the 
H.E.S.S. galactic plane survey with the flux larger than 
0.03 Crab above 1 TeV is nearly complete at present), then 
£-p,50 0.01 is implied, otherwise A'^tcV for the case (1) is 
much larger than a few. This might suggest the escape of 
high-energy particles from the SNR starts by tagc ~ 10^ yrs, 
or injection efflciency is low for old SNRs. Then, the number 
of TeV sources, NtbV, for the case (3) should be small. On 
the other hand, the case (2) may be still likely because the 
number of sources for this case can be comparable to that for 
the case (1). Hence, some TeV sources are associated with 
the GMC, but others not. In order to clarify whether the 
TeV source is associated with a GMC, the CO observation 
may be important. 
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